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Abstract 

When describing elastic deformations of a body sometimes it is 
worth to take in account elastic spatial dispersion. If spatial disper¬ 
sion is weak, as usually happens, then it can be reduced to dependence 
of thermodynamic potential on strain gradients. Such a dependence 
may be worth in association with small body size which imply large 
gradients. Besides, the inclusion of this dispersion leads to physical 
phenomena absent without it. An example of the latter is flexoelec- 
tricity. Remarkable fact is that while the derivation of differential 
equations of elastic equilibrium can be made by ordinary means in 
this case, the derivation of boundary conditions for them is less trivial 
task. This is due to the fact that strain gradients should be repre¬ 
sented in terms of second gradients of independently varied elastic 
displacements. Detailed consideration of this problem is the subject 
of this paper. 


1 Introduction 

The classical theory of elasticity is based on the thermodynamic potential 
depended on the elastic strain only. However, in some cases it is important 
to consider dependence of the thermodynamic potential on the gradients of 
elastic strain. 
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There are two main reasons why consideration of such a thermodynamic 
potential may be important. First, elastic body can be of small size, small 
sizes imply large gradients. This becomes especially important due to rapid 
development of nanotechnology. Secondly, such a dependence may lead to 
physical phenomena absent without it. An example of such phenomenon is 
the flexoelectric effect (see BE] and references therein). 

The basic problem of the theory of elasticity is to hnd static equilibrium 
of a body. Generally speaking dynamic problems are also important, but we 
restrict ourselves to the static case. Generalization of the considered theory 
to the dynamic case is absolutely straightforward. 

Even if thermodynamic potential depends on strain gradients, it remains 
quasi-local potential (it can be represented as integral of the corresponding 
density). Thus, equilibrium is described by differential equations in partial 
derivatives. This equations should be appended with boundary conditions on 
the surface of the body. Both differential equations and boundary conditions 
are determined by the conditions of thermodynamic potential minimum in 
framework of the calculus of variations. 

It is important to note that in the considered case the derivation of differ¬ 
ential equations of elastic equilibrium does not differ much from the case of 
the classical theory, it can be made by completely standard methods. On the 
contrary, the correct derivation of the boundary conditions for these equa¬ 
tions is non-trivial task. Recent studies on the subject have been done with 
regard to the description of the flexoelectric effect, but they are applicable 
in all cases when the thermodynamic potential depends on strain gradients. 

Essence of matter leading to the fact that in this case the problem of 
boundary conditions is non-trivial is as follows. Dependence of the thermo¬ 
dynamic potential of the strain gradient actually means its dependence on the 
second spatial derivatives of the elastic displacements Ui. When calculating 
the variation of the thermodynamic potential one should express it in terms 
of the independent variate. As usual, this is done using integration by parts. 
However, in the considered case surface integrals, which arises due to such 
an integration, contained 5ui^j not just 5ui. Hereinafter index separated by 
comma denotes spatial derivative in respective coordinate. Sometimes such 
derivative is also denoted as di. Values 5ui^j are not independent variates. 
Thus, additional analysis and mathematical constructions are required. 

Apparently a hrst solution to this problem has been proposed by Toupin 
[3 SI many years ago. This solution is based on the representation of the 
gradient di as a sum of “normal gradient” UiD and “tangential gradient” , 
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where rzj is unit vector normal to the surface. Operator D has been defined 
as riidi. Hereinafter the Einstein summation rule over repeated indices is 
implied. Note that formal dehnition of Di has not been proposed. Certainly, 
if Di is applied to a function, dehned in whole 3-dimensional space (at least 
in the body bulk), then Di can be well-defined as difference between di and 
UiD , no problem appears in this way. However, Toupin’s boundary condi¬ 
tions contains not only such mathematical objects, but also Bij = DiUj . 
Vector Hi is defined only on the surface and hence neither di nor D can not 
be applied to it. Thus, mathematical sense of 3-dimensional tensor Bij has 
been remained unclear. Of course, this construction can be provided with 
some intuitive sense in Cartesian coordinate system, if a couple of axes are 
tangential to the surface at given point. But anyway it is not flawless mathe¬ 
matical dehnition, because tangential plane has the only point common with 
the curved surface. Formally, vector u, is not dehned in some, at least in- 
hnitesimal, region of this tangential plane. Naturally, the derivative can be 
dehned only if diherentiated expression is dehned in some region. 

In [21, m tensor Bij has been referred as “second fundamental form” of 
the body surface. Note that second fundamental form of the surface, at least 
within conventional dehnitions, is dehned in 2-dimensional space tangential 
to the surface, rather then in whole 3-dimensional space. Thus, relation 
between second fundamental form and 3-dimensional tensor Bij remains un¬ 
clear. At least additional geometric constructions and dehnitions are required 
here to make it meaningful. 

In subsequent years, Toupin’s construction has been used to derive bound¬ 
ary conditions in many papers (see for instance). However, mathe¬ 

matically hawless dehnition of DiUj has not been given. Moreover, in these 
papers it was not even mentioned the relation between this tensor and the 
second fundamental form of the boundary surface. 

In the framework of continual theory of hexoelectricity an attempt to 
derive the mechanical boundary conditions was made also in the paper [8]. 
The authors of this paper quite correctly pointed out that the mechanical 
boundary conditions can be obtained from the surface integral arising when 
the thermodynamic potential is varied in the elastic displacements. It was 
also correctly noticed that since the thermodynamic potential depends on 
strain gradients in this case, the differential equations are of the fourth order 
in contrary to the standard theory of elasticity, where they are of the second 
order. Therefore, the number of the boundary conditions should be grater 
then such a number in the classical theory of elasticity. 
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Nevertheless mechanical boundary conditions written in [5] can not be 
considered as correct ones. Authors do not give detailed derivation of them, 
so it is not possible to specify where a mistake done. However, erroneousness 
of these boundary conditions is clear already from the fact that they do 
not contain the vector normal to the boundary surface. This means, among 
others, that it is impossible to passage to the limit of the classical theory of 
elasticity where boundary conditions contain such a vector. 

In addition it is obvious that the authors of [8] did not take in account 
the curvature of the boundary surface. The surface of the hnite-size body 
can not be flat elsewhere, further it is shown that in the considered case the 
surface curvature plays an important role. 

Number of articles [aiiniiii] have been published also where for the case 
of presence of flexoelectricity mechanical boundary conditions on the free 
surface of a body were declared in the form aijUj = 0 which is analogous to 
the form of the boundary conditions in the classical theory of elasticity. Here 
CTjj is tensor of non-classical “physical stress” which enters into differential 
equations of equilibrium in the form cr^j- = 0. Note that this boundary 
conditions were simple declared without any derivation. Note also that, as 
mentioned above, number of boundary conditions should be greater at least. 

Mathematically strict way to derive the mechanical boundary conditions 
in the presence of flexoelectricity was independently proposed in the paper 
[12]. In this paper, in contrast earlier works, projector to the body surface, 
needed to decompose Sui^j to normal and tangential parts, was constructed 
by means of two-dimensional curvilinear coordinate system on the body sur¬ 
face. This two-dimensional coordinate system was introduced additionally to 
common three-dimensional Cartesian coordinate system in the space. Such 
mathematical construction does not lead to ill-dehned mathematical objects. 
Note that although in this paper the boundary conditions has been derived 
in framework of the theory of flexoelectricity, in fact they can be applied to 
all cases where thermodynamic potential depends on strain gradients. 

Latter yet another approach to the derivation of such boundary conditions 
was proposed in the paper [T3|. This approach is based on a mathematical 
tool which is almost identical to the one used in general theory of relativity. 
The only difference is that, since the space is flat, the Riemann tensor is 
equal to zero and covariant derivatives commute. 

Note that for solution of certain problems the curvilinear coordinate^ x" 

^Note that index distinguishing curvilinear coordinates conventionally written as su- 
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are more convenient then Cartesian ones, especially if one choose a coordinate 
system, in which the equation of the body surface is , where Xg is 

some constant. In such coordinate system transformation of surface integrals 
needed for derivation of discussed boundary conditions is performed directly 
without any mathematical tricks. That is exactly what was used in [13]. 

Detailed derivation of the boundary conditions in curvilinear coordinates 
has not been given in [T3], there was only a brief note about the method 
of derivation. It was also omitted the proof that this boundary conditions 
are completely equivalent to the boundary conditions derived in [12]. This 
equivalence was only mentioned without any evidence. Consideration was 
too brief in [T3] . 

Thus, detailed description of the problem of mechanical boundary condi¬ 
tions for a case where the thermodynamic potential depends on strain gra¬ 
dients is desirable. It is presented here and in somewhat more general form 
then in dans]. 

2 General description of the mechanical equi¬ 
librium of the elastic body 

In the most general case, the mechanical equilibrium of an elastic body can 
be found in framework of calculus of variations by minimization of thermo¬ 
dynamic potential H. The minimization should be performed with respect 
to variations of the elastic displacements Ui. 

Thermodynamic potential also depends on variables other than Uj. Par¬ 
ticularly, in the theory of flexoelectricity it also depends on electric polariza¬ 
tion at least. However, here we consider only the mechanical aspects of the 
theory, so that electric polarization and all other variables, which differ from 
Ui , will be considered as given. 

In quasi-local case the thermodynamic potential of the body is the sum 
of bulk part Hb and surface part Hs- The bulk part is volume integral of 
the bulk density of the thermodynamic potential: 

HB = JnBdV. ( 1 ) 

Analogously the surface part is an integral over the body surface of the 
perscript instead of subscript. 
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surface density: 

Hs = jnsdS. (2) 

Note that the decomposition of the thermodynamic potential to volume 
and surface parts can be non-unique. If "Hb contains a term which is di¬ 
vergence then the volume integral of this term can be converted to a surface 
integral. The reverse transformation is possible also if "Hb contains a suitable 
term. 

As indicated above to hnd the body equilibrium it is necessary to calculate 
the variation of the thermodynamic potential and equate it to zero. The 
result of this calculation depends essentially on the type of the functional 
dependence of "Hb and l-is on ^i. 

In the classical theory of elasticity "Hb depends only on the strain tensor 
Uij which is the symmetric part of the distortion tensor: 

1 

Uij T Ujj') . (3) 


In the case considered here "Hb depends on Uij and also depends on Uij^k- 
Note that in the general reasoning it is quite possible to assume that "Hb 
depends on Uij and Ui j,a: directly. Real dependence only on Uij and Uij^k is 
automatically ensured due to the specific form of this dependence in this way, 
in fact it is ensured by symmetry properties of the material tensors. Keeping 
this in mind further we use unsymmetrical distortion tensor Uij . 

As for Us in the classical case it depends on Ui only, particularly this 
dependence describes interaction with the external forces acting on the body 
surface. In the case considered here one should take into account also the 
dependence of Tis on Uij. Dependence Hs on Uij is of the same order of 
smallness as the dependence "Hb on Uij^k- This is clear already from the fact 
that volume integral of "Hb being a divergence can be transformed to surface 
integral, order of derivatives is reduced at this transformation. 

For such "Hb and Tis the variation 6H can be represented in terms of 
following quantities: 

^ dy-B 


0 


ijk — 


&Hb 

dUij^k 


(5) 


Fi = 


dUs 

dui 


( 6 ) 
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( 7 ) 


dHs 


du. 




In this notation following equations can be written immediately: 


J' TijduijdV j ^ijkduij^^dV , 


( 8 ) 


6Hs = j FiduidS + ^ J:ij6uijdS . (9) 

Variations 6ui are independent but 6uij and Suij^k are not. So that 
integrals in ([8]) should be transformed using integration by parts. This yields 

dddj^ J' (T{j j ^ijk,j,k} duidV 



where rii is unit vector normal to the surface of the body. So that surface 
terms appear which are analogous to the terms in right hand side of (l9ll . 

Variations 6ui are absolutely arbitrary. Particularly 6ui and 6uij can 
be vanishing on the surface. In such a case all surface integrals disappear, 
remaining volume integral yields differential equations of equilibrium; 

~ ^ijk,j,k = 0 . ( 11 ) 

Note that if one dehnes a “physical stress” cr^ as follows 

Tij Qijk,k ( 72 ) 

then flTT]) can be rewritten in the manner analogous to the case of the classical 
theory of elasticity: 

o-ijj = 0 . (13) 

Boundary conditions for equations flT^ should be derived from the nul- 
lihcation of sum of surface integrals for 6ui and 6uij not vanishing on the 
surface. So that boundary conditions are contained in the integral equation 

j) {(TijUj + Fi) duidS + j) {Eij + QijkUk) duijdS = 0 , (14) 

where we use flT^ to rewrite the first integral. 
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In the classical theory of elasticity 0^^ = 0, Tnj = 0 so that there is no 
second integral in fll4p which contains 6uij . Therefore, in this theory there 
is no problem to derive the bonndary conditions from flTT|) at ones. One need 
only eqnate to zero the expression in parentheses in the hrst integral. 

In the case when the thermodynamic potential depends on Uij^k the sit¬ 
uation becomes more complicated. Matter is that 6ui and 6uij can not be 
considered as completely independent. Therefore, to obtain the boundary 
conditions one should transform the second integral in left hand side of (fTT|) . 
This transformation is non-trivial, in order to make it one needs a special 
mathematical tool which is considered in the following sections. Only after 
this transformation mechanical boundary conditions can be obtained. 


3 Additional surface coordinate system and 
boundary conditions in Cartesian tensor 
components 

As stated above, to derive the boundary conditions one must transform the 
second integral in the left hand side of ffTTl) . One method of this transforma¬ 
tion based on the use of additional surface coordinate system was proposed 
in [12]. In this section we consider this method with some generalization 
relatively to the paper [12]. This generalization is quite direct but it can be 
useful sometimes. 

Besides the usual Cartesian coordinate system in the space, let be an 
additional two-dimensional coordinate system x" on the surface of the body, 
a = 1,2 here. Generally the closed surface of a body is curved. Certainly 
body surface can be flat everywhere except the sharp edges where flat parts 
of surface intersect. Such a situation can be considered as the limit of the 
body with smoothed edges for which the curvature of the surface is present 
only on this smoothed edges. It is what further implied. 

In this way a two-dimensional non-Euclidean geometry appears on the 
surface of a body. This is why the surface coordinate system should be 
described as curvilinear and one should distinguish corresponding tensor in¬ 
dices, denoted by Greek letters, as superscripts and subscripts. Indices cor¬ 
responding three-dimensional Cartesian coordinates are denoted by Latin 
letters, they are always subscripts. 

If only two-dimensional tangential to the surface tensors were necessary 


then the well-known mathematics of differential Riemannian geometry would 
be sufficient. However, Eq. fll4p contains three-dimensional tensors on the 
surface rather then two-dimensional tensors. This is why the standard tech¬ 
nique of Riemannian geometry should be modihed here. Such a modihcation 
is quite natural if the frame formalism is used. 

Let us introduce the frame vectors 

eia = dri/dx"^ = , (15) 

where are Cartesian coordinates of a point on the surface. Consequently, 
if point moves along the surface then 

(dr)^ = eiaCipdx'^dx^ . (16) 

So that metric tensor of the surface coordinate system can be represented 
as follows 

9aP ^ia^ip • 

Since generally the coordinates are curvilinear, the basis set Cia is 
neither orthogonal nor normalized. Therefore it is necessary to introduce a 
mutual basis set ef, which is dehned by relations 

= d^p , (18) 

where 5^ is Kronecker delta in Greek indices. By means of the reverse metric 
tensor mutual reference vectors ef can be represented in terms of eia- 

e“ = . (19) 

Thus, by means of gap and g'^^ curvilinear indices can be raised and 
lowered in the standard manner. Surface integrals can be represented in 
the form of integration over a;" using obvious relation: dS = y/gd^x, where 
d‘^x = dx^dx’^i g = detgap- 

The key observation needed further is that e^Cja is the three-dimensional 
projector onto the plane tangential to the surface of the body. Indeed, 
efejaejCkp = efd^Ckp = efcka- Moreover any vector Vj convolved with efcja 
yields a linear combination of vectors e" which are tangential to the body 
surface. These two facts prove our assertion. It is also obvious that niUj is 
the projector onto normal of the surface. Thus, the unit operator (Kronecker 
delta in Latin indices) can be represented as follows: 

Sij TiiTlj -|- 
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Eq. fEUD allows us to convert Sutj in such manner: Suij = dui^inirij + 
Sui^iCiae'^. From Eq. (ITB]) it follows that Sui^iCia = where (.. .)_o, 

means the partial derivative in Therefore, 


(.. .)6uijdS = ji.-- ■)Sui^ininjy/g(fx + ^(- • .)ej{6ui)^ay/gd‘^x . (21) 

The second term in the right hand side of this equation can be integrated by 
parts in straightforward manner. It yields 


.. .)5uijdS = f{...)6ui^ininjdS- f [(. ..)e“v^j 5uig ^^“^dS. (22) 


Identity fl2^ allows to transform the surface integral contained not in¬ 
dependent variations 6uij to the surface integrals contained 6ui and dui^iui 
which are independent. Indeed, 5ui are independent initially. As for dui^iui 
it is nothing but the gradient in the direction normal to the surface. Setting 
6ui on the surface one can continue them arbitrary along the normal to the 
surface. So that near surface Sut^ini are arbitrary independent values. 

Thus, by applying (|2^ to flTT)) and taking into account the independence 
and arbitrariness of Sui and , we obtain the boundary conditions in 

the form 

aijUj + Fi- ^ - (QijknkeJg^^'^) ^ = 0 . (24) 

Naturally, these equations are valid only on the body surface. They are 
boundary conditions for differential equations of elastic equilibrium flT^ . 

Using the Leibniz rule, identity Qijk,a = Qijk,Fi,ay Eq. (IT^ and Eq. (1201) 
the last term in left hand side of fl2T)) can be represented in the form used in 
the paper [12] : 

(Qijkrike^g^^'^) ^ = Qijkjnk - &ijk,ininjnk + Qijk'ljk , (25) 


where 

7ifc = ^ . (26) 

Generally an analogous transformation of rnay be 

meaningless because Sjj may be dehned only on the surface rather than 
space. However, there are special case when Sjj can be represented as con¬ 
volution: 

Eij = Sijiui , (27) 


10 




where tensor Siji is well-defined in the bnlk of the body. In particnlar, snch a 
possibility appears when the snrface contribntion to the thermodynamic po¬ 
tential arises dne to the transformation of volnme contribntion by integration 
by parts. For snch cases can be transformed absolntely 

similarly as (QijkUkejg^^^^ 

Eqnations fl25]l . fl21]) with the eqnalities fl2^ . are nothing bnt some 
generalization of the bonndary conditions obtained in [T2]. Generalization is 
as follows. First, we have consider the general case of the thermodynamic 
potential depended on the strain gradient (see previons section) rather then 
a special thermodynamic potential nsed in |T2] . Second, we have taken into 
acconnt the snrface contribntion to this potential, it has led to additional 
therm in 0231) and additional terms F) and — in 

([21]). Except for these generalizations bonndary conditions derived above 
coincide with ones derived in jl 2 j . 

Note that bonndary conditions fl23D . fl 2 ll) do satisfy the condition that in 
the limit case they should correspond the classical theory of elasticity. Indeed 
in the classical limit Sjj = 0, 0*^^ = 0. So that 0231) becomes an identity 
while 0241) takes the classical form (Tijfij + Fi = 0. 

Remarkable fact is that curvature of the body surface plays an important 
role in the mechanical boundary conditions derived above. To see it let us 
rewrite fl2l|) by means fl25l) omitting surface contribution contained It 
turns out: 

^ij^j T ^ijk,j^k T ^ijk^jk 0 ■ 

Tensor 'yjk is a purely geometric construction (see also Appendix]^. It does 
not depend on the elastic displacements or somewhat else of non-geometrical 
nature, actually it characterizes only the geometry of the body surface. Note 
also that this tensor is invariant with respect to the change of the surface 
coordinate system. It is clear from the fact that it is a scalar convolution in 
Greek indices. If some part of a body surface is fiat then suitable Gartesian 
coordinates can be chosen as surface coordinates on this part. But in such 
surface coordinates all values under derivative in (12^ are constants. So that 
-yjk = 0 on the fiat parts of the body surface. Otherwise -yjk 7 ^ 0 on the curved 
parts of the body surface, it is absolutely obvious. Thus, the last term in 
left hand side of fl28|) appears or disappears depending on whether the body 
surface is curved or fiat. A similar conclusion can be drawn regarding the 
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surface term g , one should only transform it by means of 

Leibniz rule. 


4 Curvilinear coordinates and covariant for¬ 
malism 

In the previous section the mechanical boundary conditions were derived 
in terms of tensor components related to the three-dimensional Cartesian 
coordinate system. To do this one has introduce an additional curvilinear 
two-dimensional coordinate system on the surface of the body. 

One can go the other way and just introduce curvilinear coordinate system 
in the whole space. If it is required that the equation of the body surface 
in this curvilinear coordinate system has a simple form = a;|, where x| is 
some constant, then the derivation of the boundary conditions can be done in 
a more direct way. It is essential that such a coordinate system additionally 
is most convenient for the practical solution of differential equations of elastic 
equilibrium. 

When curvilinear coordinates are used the space remains flat. However, 
to work with tensor objects at these case is most convenient in framework 
of the mathematical tool of Riemann geometry well known from the general 
theory of relativity. The fact that space is still flat only leads to that the 
Riemann tensor is identically zero. 

Although the above-mentioned mathematical formalism (so called covari¬ 
ant formalism) is well known in theoretical physics, it is not widely used in 
material science. This is why we describe its main points in this section for 
convenience. It is also needed for references. 

It is convenient to assume that besides curvilinear coordinate system there 
is Cartesian coordinate system also. So that we can introduce frame vectors 
Cia by the same equation (IT^ from the previous section. The only difference 
is that now Greek indices run three values 1, 2 and 3 and now means 
Cartesian coordinates of any point of space. Taking into account this com¬ 
ment the equations flTT)) . flTK]) and (IT^ from previous section also remain 
valid. Contrariwise equation fl20|) now is changed into 

efeja = Sij . (29) 

The curvilinear contravariant tensor components are dehned by Cartesian 
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tensor components as follows: 


Aai...an — A . . 

^ ~ 


(30) 


The cnrvilinear covariant tensor components are defined analogonsly: 




(31) 


By analogy it is clear how are defined mixed tensor components. By means 
(12^ it is obvions how cnrvilinear tensor components can be transformed back 
to Cartesian tensor components. Say for a tensor in covariant components 
this transformations looks as follows: 


(32) 

Tensor convolntions are transformed by means of fflSj) . For example 

. (33) 

Thns it is obtained antomatically that only npper and lower Greek indices 
can be convolved. 

Gradients of tensors are transformed in a special way. Since the frame 
vectors are not constant they can not be moved throngh the derivative in 
straightforward manner, additional terms appear dnring this motion. For 
example gradient of vector is transformed as follows: 

where 

Fq,^ G G • (2^) 

The expression in the right hand side of fIMl) is called the covariant derivative 
and is denoted by index separated by a semicolon: 


^a;/3 — ^a,/3 ~ ■ (36) 

Set of valnes F^)^ is called the Ghristoffel symbols. They can be expressed 
not only in terms of frame vectors (l3^ bnt also in terms of metric tensor 



g'^^{gsa,l3 + g5^,a 


gal3,5) ■ 


(37) 
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It is easy to see that this Christoffel symbols are symmetric in the lower 
indices. There is also useful equation 

h5 = |^ = (lnV9),« (38) 

which can be obtained if one differentiate In = Sp In gap , g = detgap here¬ 
after. 

The covariant derivative of contravariant vector should be dehned in an¬ 
other manner. The reason is that frame vector e“ should be moved through 
the derivative rather then Cjo. Note that is a constant (see equation 

01 ), so that it partial derivative is zero. Therefore and 

from this identity and simple algebra it turns out: 

Mjeyu, = . (39) 

From the above it is quite clear how to determine the covariant derivative 
of higher rank tensors . The terms containing the Christoffel symbols (so- 
called T-terms) appear due to frame vectors motion through the derivative. If 
the tensor is of the highest rank then several frame vectors should be moved 
through derivative. Obviously several T-terms should appear in this way. 
For example covariant derivative of two-rank tensor with lower indices is 

^a/3;7 = ^0/3,7 “ “ ^aS^^p^ ■ (40) 

What happens in the other cases is clear by analogy. Covariant derivative of 
a scalar is determined as usual partial derivative. 

It is easy to show that the covariant derivative of the product obeys the 
Leibniz rule 

(^....B...);o = A B -a ■ (41) 

It can be also proved in straightforward manner that the covariant derivative 
of frame vectors is identically zero. Since metric tensor can be represented as 
the convolution (in Cartesian indices) of frame vectors, it covariant derivative 
is identically zero also. It is also worth that if some tensor is constant in 
Cartesian components (say it is a material tensor of homogeneous media) 
then in curvilinear components it covariant (but not usual partial) derivatives 
are zero. 

In the general case of curved space covariant derivatives do not commute, 
their commutator is proportional to the Riemann curvature tensor. However, 
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here we are interested only in the case of curvilinear coordinates in flat space. 
In the flat space Riemann tensor is identically zero. Therefore, the covariant 
derivatives commute. 

To complete this brief overview of covariant tensor analysis it remains to 
say few words about integration. As for the volume integration everything is 
quite simple. It is easy to see that the Jacobian of the transformation from 
Cartesian coordinates to curvilinear ones is the determinant of Cia ■ However, 
it is more convenient and generally accepted to express this Jacobian in terms 
of metric tensor. It is clear from (fT7|) that the square of the Jacobian equal 
to g. So that 

dV = y/gd^x , (42) 

where d^x = dx^dx'^dx^. This formula is enough for volume integration, at 
least in the part that we need here. 

As for the integration over the surface, the situation is somewhat more 
complicated. Here we restrict ourselves to the case when the equation of the 
surface has a simple form x^ = const. The more general case we do not need. 

If the equation of the surface has the simple form mentioned above then 
the integration over the surface reduces to integration over x^ and x^. Let us 
introduce additional notation: Greek indices in parentheses have the values 
1 and 2 only. In this notation it turns out 


dS = y/g d^x , 


(43) 


where d?x = dx^dx "^, g is determinant of 2 x 2 matrix g(a){i 3 ) ■ Note that 
according to the known representation of the inverse matrix by algebraic 
adjuncts 

g = • (44) 

Surface integrals often contain an unit vector normal to the surface. So 
that we also need the expression for such a vector. It is clear geometrically 
that in the considered special coordinate system the vectors en and 6*2 are 
tangential to the surface. Therefore in accordance with (ITSll the vector ef is 
normal to the surface. By making the normalization of this vector, we obtain 
the unit vector normal to the surface in Cartesian components: 


rii = 






(45) 
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Conversion of this vector to the curvilinear components by means of general 
equations yields the following: 


Tl/y - TljCi, 




(46) 


a3 

n- = g-% = ^. ( 47 ) 

The above brief consideration of the covariant formalism is sufficient for 
the purposes of subsequent exposition. For more details one should refer to 
any textbook on general relativity. 


5 Boundary conditions within covariant for¬ 
malism 

In framework of mathematics briefly described in the previous section equa¬ 
tions flT^ - flTTll take the form 

= (49) 

i bUa\fz‘^x + f — 0. (50) 

One should stress that equation of the body surface of the form = 

const is implied here, otherwise the last equation is not valid. Should be 
noted also that the differentiation of scalar in covariant tensor components 
yields contravariant tensor components. So that 


j.a0 _ dy-B 

(51) 

0a^7 _ dT~LB 

^^ck;/3;7 

(52) 

po. _ dUs 
dUa 

(53) 

dlic)L-Q 

(54) 
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Transformation of the second integral in left hand side of fISOp can be pro¬ 
vided quite directly. One should represent 6 ua-i 3 in terms of partial derivative 
and T-term 

and separate the case of /3 = 3 in the terms contained 5ua,p as follows: 

(.. = (.. + (.. . (56) 

Remember that index in parentheses means the index which values are 1 and 
2 only. Integral contained 5ua,{p) can be integrated by parts without any 
tricks. In this way it turns out: 

j> + j) -I- 6ua,3y/g(fx = 0 , (57) 


where 

L“ = + F“ - - (E“*'’>g‘''“) 

(58) 

Note that since (jH]) and (H6P L" can be rewritten using g rather then g in 
last term: 


L“ = + F“ - - (E“l'’>g'''^) g“‘''^- 

(59) 

From flFzp the boundary conditions are obvious and have the following 
form 

S“3 + ©“37^^ = 0 , (60) 

L“ = 0. (61) 

Note that since rvj y ( see fl46p l fl6Up can be rewritten in the form 

+ ©“^T'n^n.^ = 0 (62) 


which corresponds to fl25p . As regards (IMP with (I55p or (15^ it looks different 
then (l2Tp . Nevertheless it can be proved that (IMP is equivalent to fl2Tp . The 
proof follows below. 
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Before the description of the proof itself we should discuss some geomet¬ 
rical matter. The curvilinear coordinate system used in this section is such 
that equation of body surface is = const. Thus, hrst two coordinates 
of this coordinate system may be well used as surface coordinates system 
used in section |3l Nevertheless, in this case one one should distinguish some 
objects depending on whether are used as hrst two coordinates of spatial 
coordinate system or as the surface coordinate system in manner of section Ei 
For such a distinction here we add additional tilde to objects related to the 
surface coordinate system. 

Not that ej(Q) = ei(^a) but This is consequence of the fact that 

g{a){p) ^ g{a){p)^ gy ] 2 ieans of simple geometrical reasoning one can see that 



„(a)3 

I _ 

QQ ^2 ' 
gi.i 


(63) 


Note also that 

g = g. (64) 

Now one can turn to the proof itself. Let us denote last two terms in (|2TD 
as . By means of fl63l) and fIMl) we can rewrite it in terms of 3-dimensional 
coordinates: 


-h 


& 




g-'/' + ( S 


(a)3 N 

.y_p3„l/2 


g33 ^3 


e,g 


g 


- 1 / 2 . 








(65) 


g 


- 1/2 


-(«) 


Second and fourth terms in right hand side of this equation can be combined 
into one term 




„ « 3 , ■ 

I _£ 1/2 

^33 ® 


r-1/2 


( 66 ) 




Since ef ~ n* the expression in parentheses is zero due to ([23]). Moreover, 
since it is zero everywhere on the surface (i.e. for every x^^), its derivative 
in is zero also. So that fl6^ is zero i.e. second and fourth terms in fl6^ 
cancel each other out completely. Thus, 

(^i = - (Sy ef ^gi/2^ g-^/2 _ g’^/^ _ 
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Further one can easily calculate contravariant representation of vector It 
turns out: 


(S.j + - 

g-i/^ _ g-1/., 


( 68 ) 


In the hrst term of right hand side of this equation range of the index /3 can 
be expanded to complete set 1,2,3. This is consequence of fl2^ and the fact 
that ef ~ rij. Further one should note that 


— p" /), p^ — p“ pirp'^p^ — — r“ P'^P^ 




(69) 


With this identity equation fl68|l takes the following form 




«(/3)o-V2 


m 


.-1/2 _ (0«(/3)7^^gl/2 


r-1/2 


(70) 


,(/3) 


Note that range of index [3 in hrst term of right hand side of this equation 
can be extended to complete set 1,2,3. This is a consequence of the equation 
fl60|) . With this extension further elementary transformations yield exactly 
dni]). Proof is completed. 

Also it is worth to transform obtained boundary conditions (l60|) and fl6T|) 
to the forms used in papers In this papers only case of F" = 0 and 

— g considered, the same is assumed while the transformation. So that 
by means of (IT^ equation (1^ can be rewritten as 


0«33 _ Q _ 


(71) 


By means of Leibniz rule, equation (138|) and representation of partial deriva¬ 
tive in terms of covariant derivative last term in right hand side of fl59p can 
be transformed as follows 






(72) 


With this identity (16T]) takes the form 




(0“37^i/2) 


7-1/2 


n.y = 0. 


(73) 
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This coincides with (12) from [T3]. By means and fl7ip it can be also 
rewritten in following form 

^a3_0a/33 ^ 0a(/35)p3^^^ ^ 0a33^^ _ q ^ ^ 74 ^ 

where one pair of parentheses enclose several indexes mean that these in¬ 
dices are not equal 3 simultaneously. Boundary conditions fITT]) and fl71|) are 
exactly the ones used in 

6 Conclusion 

It is described above how one can generalize classical theory of elastic con¬ 
tinuum to a case when thermodynamic potential depends on the strain and 
the strain gradients. In this case, in fact, one gets such a type of classical 
(non-quantum) held theory were function which determinates held equations 
(“Lagrangian”) depends on the hrst and second derivatives of the initial 
held, particularly the elastic displacement held Ui. Certainly in the static 
case “Lagrangian” is nothing but the thermodynamic potential. 

At the same time it should be noted that in essence we did not use the fact 
that the held Ui is a held of elastic displacements. All this is also true in cases 
where the vector held Ui has another physical nature if only “Lagrangian” 
has the form that we has expected. It is not a signihcant problem to add 
also a “mass term” (in terms of the quantum held theory) aijUiUj where aij 
is some constant tensor. It will not change much, only diherential equations 
will be slightly (and obviously) modihed while the boundary conditions will 
remain the same. Such opportunities should be kept in mind, but practically 
considered theory is needed to describe namely the elastic helds, and mostly 
in relation to the hexoelectric ehect. 

It is shown in this paper that in theories of this type, particularly in 
generalized theory of elasticity, diherential equations can be derived in quite 
conventional manner. The only diherence with respect to ordinal case is 
that while varying one should integrate by parts twice. As for the derivation 
of the boundary conditions for these diherential equations, the situation is 
more complicated. We have considered two ways by which such derivation 
can be done. The hrst way uses a mathematical trick based on the use of 
additional surface coordinate system. The second way is more direct, but it 
involves the use of a covariant tensor analysis and a curvilinear coordinate 
system of a special class. In the coordinate systems of this class equation of 
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boundary surface should have the form where x| is some constant. 

For a complex surface shape it may be difficult to hud a single coordinate 
system of this type, but in this case it is possible to divide the surface to 
such parts, for which the hnding of the coordinate system of the desired type 
is not very difficult. It is, however, necessary to ensure that these parts are 
sewed sufficiently smoothly. 

It should be noted also that in the considered case the differential equa¬ 
tions and especially boundary conditions are very complicated. With regard 
to the calculation of the flexoelectric deformations corresponding boundary 
problem was solved for a homogeneously polarized ball of an isotropic ma¬ 
terial [13]. Even in this simple geometry one has to use a representation 
of the radial functions in the form of power series. To solve more complex 
boundary problems even more difficult. But for the calculation of specih- 
cally flexoelectric deformations it was proposed an approximate method [13] 
which radically simplihes the problem. Corresponding theory is considered 
in detail in another paper [15] . 
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Appendix A Relation between tensor 7 ^j and 

Gauss second fundamental form 


By means of Leibniz rule equation fl26l) . which defines the tensor 7 ^ , can be 
rewritten as follows: 


lij 


n 




ef + n^e" 


i,a + '^j^i 


(75) 


It can been shown (see Section 01) that 


9,oi 99^^^ 913^,a 


(76) 


Representing metric tensor and reverse metric tensor in terms of surface 
frame vectors, we hnd: 

9,a ‘^9^1 

This equation can be essentially simplified. By means of the relation = 
6ki — riiUk , it turns out 

= ‘^9 (efefc/3,« - e^ek^^aniUk) ■ (78) 

Any surface frame vector is orthogonal to vector ni normal to the surface. 
Thus, second term in parentheses is identically zero and it turns out: 

9,a = 2gelekf5,a ■ (79) 

Thus, fITHll can be rewritten as follows: 

lij = rij^aCi + njel^ + nje^eiekg,a ■ (80) 

Second term in right hand side of (1801) can be transformed in such manner 
that its appears the term cancelling out third term. First, we transform it 
as follows: 

= njSikel^ = Uj-ef Ckpe^^ + ritrijUkel^ . (81) 

From the fact that ekjse-k^ = does not depend on x" , and obvious relation 
efc/ 3 ,o = Cka^p , it follows that ek^e'^a = ~^ka,i 3 e^ ■ Thus, it turns out: 

= riiTijUkel^^ - Ujefekfi^aei ■ (82) 
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Last term in right hand side of this equation obviously cancels out third term 
in right hand side of flSU]) . Thus, 

lij = + niTijUkc'^^^ . (83) 

First term in right hand side of fl83ll can be related to Gauss second 
fundamental form ba/s , which is dehned as follows: 

bal3 = nkekl3,a ■ (84) 

Since nkCkp = 0 , Uk^a^k/s = —nk^kp^a ■ Convolving this equation with and 
taking into account that Ckpe^j = 6jk — fijUk , we hnd: 

flj^a ^k,a^j^k ^k(^kp,a^j • (85) 

From nkUk = 1, it follows that Uk^a'^k = 0. Thus, 

7T.J Q '^k^kp,a^j , (86) 

In the theory of surfaces this equation is known as Weingarten equation. By 
means of it we hnd: 

= -Bij , (87) 

where 

Bij = bape^e^ . ( 88 ) 

It should be emphasized that 3-dimensional tensor Bij is not a tensor of 

second fundamental form bap itself, although it is related to b^p . 

Second term in right hand side of fl83p can be also represented in terms 
of Bij . Obviously Bn = bapg^^ ■ Next, 

niTijUkek^^ = ninjnkg "^^-7 ninjnkg°‘^ekp,a ■ (89) 

First term in right hand side of this equation is zero identically, because 
rifcef = 0. It remains only the term, which contains UkCkp^a = b^p . This 
eventually yields the trace of the tensor B^j . Finally it turns out: 

Fp ^i'^jBkk * (^6) 
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